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We show that Wronskians between properly chosen linearly independent solutions of the 
Schrodinger equation greatly facilitate the study of quantum scattering in one dimension. They 
enable one to obtain the necessary relationships between the coefficients that determine the asymp- 
totic behavior of the wavefunction. As illustrative examples we calculate the transmission probability 
for the penetration of a Gaussian barrier and the scattering resonances for a Gaussian well. 
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I. INTRODUCTION 



Barrier penetration or tunnel effect is one of the most striking predictions of quantum mechanics and most textbooks 
discuss the few available exactly solvable one-dimensional models^—. In this journal there has been great interest in 
the subject^— as well as in one-dimensional potential scattering in genera l 7 ' 8 ' 17 " — . Such interest has been focussed 
not only on the analytical properties of potential scatterin g 17 i 19 i 20 but also on the numerical calculation of transmission 
probabilities and other quantities that describe barrier penetration and potential scatterin g 9 ' 11 ' 12 ] 14-16 ] 21 . 

The purpose of this paper is to show that a straightforward application of Wronskians, which are well known in 
the study of ordinary linear differential equations^ 2 ^, greatly facilitates the discussion of the analytical properties of 
potential scattering. They are also suitable for the numerical calculation of the transmission probability and any other 
quantity of interest in potential scattering. Present discussion is motivated by an earlier application of Wronskians 
to the analysis of one-dimensional models for resonance tunneling reactions 2 ^. 

In Sec. [II] we apply Wronskians to the Schrodinger equation for potential scattering in one dimension and derive 
matrix equations connecting the coefficients of the asymptotic forms of the wavefunction left and right of the scattering 
center. In Sec. IIIII we specialize in a general short-range interaction and derive equations for the transmission 
probability. As illustrative nontrivial examples we consider a Gassian barrier and a Gaussian well. By means of an 
exactly-solvable problem we discuss the origin of the scattering resonances in potential wells. In Sec. II VI we summarize 
the main results and draw conclusions. In order to make this paper sufficiently self-contained we collect some well 
known mathematical properties of the Wronskians in an Appendix. 



II. WRONSKIANS AND THE SCHRODINGER EQUATION 

The time-independent Schrodinger equation for a particle of mass m that moves in one dimension (— oo < X < oo) 
under the effect of a potential V(X) is 

-^"{X) + V{X)^{X) = E^{X) (1) 
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If we define the dimensionless coordinate x = X/L 7 where L is an appropriate length scale, then we obtain the 
dimensionlcss eigenvalue equation 

+ v(x)ip(x) = eip(x) 

(p(x) = y/Lip{Lx), v(x) = ^-V(Lx), e = (2) 

The length unit L that renders both e and v(x) dimensionless is arbitrary and we can choose it in such a way that 
makes the Schrodinger equation simpler. We will see some examples in Sec. IIIII 

In most cases of physical interest we can write the asymptotic behavior of the dimensionlcss wavefunction (p(x) as 
follows 

Aid(x) + BiSi(x) for x -)• oo 
(p(x) ^ \ (3) 
I A 3 C 3 {x) + B 3 S 3 {x) for x ->• -oo 

and in the intermediate region — oo < x < oo we have 

<p(x) = A 2 C 2 (x) + B 2 S 2 (x) (4) 

The form of the functions Cj (x) and Sj (x) depends on the problem and throughout this paper we choose all of them 
real. On the other hand, the coefficients Aj and Bj may be complex. We discuss some examples in Sec. IIIII for the 
time being we assume that those functions satisfy the conditions (|27|) at a conveniently chosen point xq. By means 
of the Eq. (|29|) given in the Appendix we can easily obtain matrix expressions connecting the coefficients Ai and Bi 
in the three regions: 





(5) 



where 



Ri 



WiC^Sx) W(S 2 ,S 1 ) 
W(d,C 2 ) W(C U S 2 ) 



W(C 3 ,S 2 ) W(S 3 ,S 2 ) 
W(C 2 ,C 3 ) W(C 2 ,S 3 ) 

In these equations W(f, g) denotes the Wronskian of the functions f(x) and g(x)&2£ already defined in the Appendix. 
If we repeat the procedure and obtain the inverse relations we appreciate that 



Rr 1 = 



W(C U S 2 ) -W(S 2 ,S!) 
-W{Cx,Ci) W{d,Sx) 



I W(C 2 ,S 3 ) -W(S 3 ,S 2 ) . 
-W(C 2 ,C 3 ) W(C 3 ,S 2 ) 
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from which it follows that 

R*J = JRJ 1 , J = (°_ i J | (8) 

where the superscript t stands for transpose. Since the skew-symmetric matrix J satisfies J* = J -1 = —J we conclude 
that 

R'J = JR 1 (9) 

which resembles the symplcctic condition for a canonical transformation in classical mechanics^. Any matrix that 
satisfies Eq. © is said to be symplectio^i. Besides, Eq. (|7|) tells us that the determinant of every symplectic matrix 
Ri, R2 and R is unity. 

If the potential is parity-invariant (y(—x) = v(x)) then C 2 {— x) = C 2 (x), S 2 (—x) = —£2(2;) and 

W{C 2 ,Sr) = W(C 2 ,S 3 ), 

W(S 2 ,S 1 ) = W(S 3 ,S 2 ) 

W(C U C 2 ) = W(C 2 ,C 3 ) 

W(C 1 ,S 2 ) = W(C 3 ,S 2 ) (10) 

so that we need to calculate half the number of Wronskians. It is worth noting that Wronskians containing C\ and 
S\ are constant for x — > 00 and those with C 3 and S3 are constant for x — > —00. 



III. POTENTIAL SCATTERING 

We assume that 

lim v(x) = v_ 

lim v(x) = v + (11) 

x — >oo 

where v± are finite constants. If v(x) approaches those limits sufficiently fast then we know that the asymptotic 
behavior of the solution is 



<p(x) -> 



A[e lklX + B[e- tklX for x -> -00 
A' 3 e iklX + B' 3 e~ ikix for x -> 00 



fci = y/2(e-v-), k 3 = ^2{e-v+) (12) 



provided that e > max{v + , 
If we choose 



C s = cosifyx), Sj(x) = Sin f jX) , j = 1,3 (13) 



in Eq. ([3]) and compare it with Eq. (fT2]) we obtain 
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Suppose that we want to study the scattering of a particle that comes from the left (x < 0). In such a case B' 3 = 
and the transmission probability is given by 

h\A'i\ 2 { ' 

Noting that B 3 = ik 3 A 3 we can easily rewrite the transmission probability in terms of the coefficients Aj and Bj 
because 



A' 3 2hA 3 



(16) 



A[ kxAi - iBi 

In the intermediate or scattering region we write (p(x) as in Eq. (j4]) where C 2 (x) and S 2 (x) arc two 

solutions of the dimensionless Schrodinger equation that satisfy Eq. ([27)1 . In order to obtain them we 

may resort to any available numerical integration method, like, for example, Runge-Kutta3£ (see also 



http://en.wikipedia.org/wiki/Runge-Kutta_methodsl. Let y{x) be either 62(2;) or 5*2(2;). Typical numerical inte- 



gration methods yield y(x) at a set of coordinate points x n ± jh where j = 0, 1, . . . and h is the step size. They 
simultaneously provide the derivative of the function y'{x) at the same set of points so that the numerical calculation 
of the Wronskians between the intermediate solutions and the asymptotic ones is straightforward. Numerical inte- 
gration methods like Runge-Kutta are available in many commercial and free softwares so that it is unnecessary to 
write a computer program for that purpose. In the present case we resorted to the fourth-order Runge-Kutta method 
built in the computer algebra system Derive ( http: / /www.chartwellyorke.com/derive.html[ ) . The starting point of the 
integration process requires 2/(20) an d 2/(2:0) that are already known for the functions 62(2;) and S 2 (x). We propagate 
the solution left and right till the Wronskians appearing in the matrices ([6]) are constant within a given error. Then we 
calculate the coefficients A\ and B\ in terms of A3 by means of Eq. (|5|) and finally the transmission coefficient from 
equations (|15[) and (|16[) . Note that the application of the numerical integration method is straightforward because all 
the functions Cj (x) and Sj (x) are real. 

In the simplest case of a parity-invariant potential we need half the Wronskians in order to obtain the matrix R in 
equations ([5]) and ([6]) as discussed at the end of Sec. HH The transmission probability reads 



i(W(C 2 , C3) 2 + kfW(C 2 , S3) 2 ) (W{C 3 , S 2 y + k 2 W(S 3 ,S 2 ) 2 )} 



and we only have to integrate the differential equations for C 2 {x) and 52(2;) from 2:0 = to the right: Xj = jh, 
j = 0,l,...,N. 

As an illustrative example we choose the Gaussian barrier 

V(X)=V e- aX \ V o ,a>0 (18) 

If we set L = 1/ \fa then v(x) = voe~ x and the dimensionless Schrodinger equation depends on just one potential 
parameter vq = mVo/(h 2 a). Since i>_ = v + = then k\ = 2e, where e = mE/(h 2 a) is the dimensionless energy 

Throughout this paper we choose the integration step size h = 0.01 and N — 500 integration points so that the 
maximum coodinate value is 2:500 = 5. Fig. [1] shows the behaviour of the functions C 2 {x) and S 2 {x) for e = 1 and 
the Gaussian barrier with vq = 2. In Fig. [2] we appreciate that the Wronskians approach constants as |x| — > 00 and 
that x = 5 is large enough for an accurate estimation of those limits. We thus calculated the transmission coefficients 
shown in Fig. [3] for three values of vq. As expected there is tunneling for all e > and lim e _ ! . 00 T = 1. It is clear that 
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the oscillatory behavior of T(e) found by Chalk^ is due to the truncation of the Gaussian potential with the purpose 
of connecting the power-series solution for the intermediate region with the asymptotic plane waves^. We appreciate 
that the necessary truncation of the integration interval does not produce any undesirable effect on the transmission 
probabilities calculated in terms of Wronskians. In fact, the great advantage of the Wronskian method is that we 
calculate the constant asymptotic limit of each Wronskian with a given desired accuracy as shown in Fig. [2] 

In passing we point out that there is no trace of the questionable tunneling condition derived by Nandi^ for the 
Gaussia barrier. 

Another simple, nontrivial, and most interesting problem is the Gaussian well 

V(X) = ~V e- aX \ V Q ,a>0 (19) 

2 

Proceeding as in the preceding example we obtain the dimensionless potential v(x) = —voe~ x and the same expres- 
sions for vq and e. In this case it is most instructive to calculate the transmission probability T in terms of vq for fixed 
values of e in order to reveal the scattering resonances as shown in Fig. |4] Note that the maxima of the transmission 
probability T = 1 that occur at some particular values of Vq are roughly independent of the energy e. This well known 
phenomenon is better understood by means of an exactly solvable model. 

The scattering resonances appearing in Fig. |4] are similar to the ones exhibited by the exactly solvable well 

V W = ~ e^AT'^^ (20) 
that we easily transform into the dimensionless potential v(x) = — vq/ cosh(x) 2 . Both the dimensionless potential 
parameter vq and energy e have the same expressions as in the preceding examples. In this case the transmission 
probability reads^ 

, 2 



sinh (ny/2e) 1 . . . , v 

T= — — y - 2 >- , A = - (1 + x/1 + 8^) (21) 

sinh (7rV2e) + sin(7rA) 2 1 



The family of curves T vs. vq for constant e resembles the one shown in Fig. |4] Note that there is full transmission 
T = 1 when sin(7rA) = 0; that is to say, when A is an integer. In order to understand the origin of these resonances 
we pay attention to the bound states^ 

e„ = -i(A - 1 - n) 2 , n = 0, 1, . . . < A - 1 (22) 

We appreciate that full transmission takes place when one of the excited bound-state energies (n > 0) lies exactly 
at the rim of the well e n = in which case A = n + 1 (the ground-state energy e = —(A — l) 2 /2 is negative for all 
vo > 0). 

This exactly solvable problem also proved to be useful for testing the accuracy of our programs for the calculation 
of the transmission probability. 



IV. CONCLUSIONS 



In this paper we propose an alternative way of approaching quantum scattering in one dimension. We think that the 
method based on the Wronskians between linearly independent solutions to the Schrodinger equation is preferable to 
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other approaches. The relatively light effort necessary to master a few mathematical properties of the Wronskians pays 
generously when attacking the scattering problem either analytically or numerically. In this paper we focused mainly 
on the latter because we are interested in nontrivial problems that are not so widely discussed in most textbooks on 
quantum mechanics^—. 

The derivation of all the necessary scattering equations in terms of Wronskians is straightforward as well as their 
practical application by means of extremely simple computer programs. The calculation of the transmission probability 
is quite reliable if one simply checks for constant Wronskians before truncating the propagation of the solutions towards 
left and right in the numerical integration routine. Thus, the error due to a finite integration interval is simply the 
error in the lack of constant Wronskians that is easily bounded to the desired accuracy. In this way we avoid any 
spurious oscillations in the transmission probability as discussed in Sec. Mil 

In addition to all that, the Wronskians have proved to be most useful for the estimation of the complex energies 
that describe tunnel resonances as discussed in the paper on which we based present pedagogical presentation^. One 
can also apply the Wronskian method to bound states and calculate their energies by simply taking into account 
the appropriate asymptotic behavior of the wavefunction and requiring that it be square integrable. It is clear that 
the approach exhibits a wide variety of useful applications and for that reason we think that it is worth teaching in 
advanced undergraduate or graduate courses on quantum mechanics. 



In order to make this paper sufficiently self-contained in this appendix we outline some well known results about 
the Wronskians that are useful for the study of ordinary differential equations in general^ and also for the treatment 
of the Schrodinger equation in particular—^. To this end, we consider the ordinary second-order differential equation 



V. 



APPENDIX 



L(y) = y"(x) + Q(x)y{x) = 



(23) 



If j/i and y 2 are two linearly independent solutions to this equation then we have 



ViL{V2) - yiL{y\) = -i-W(y 1 ,y 2 ) = 



(24) 



where 



W(yi,y 2 ) = yiy' 2 - ViV\ 



(25) 



is the Wronskian (or Wronskian determinant^). Two obvious properties are: 



W(f,g) 



W(gJ), W(f,f) = 



(26) 



By linear combination of y\{x) and y 2 (x) we easily obtain two new solutions C[x) and S(x) satisfying 



C(x ) = S\x ) = 1, C'{x ) = S{x ) = 



(27) 



at a given point xq so that W(C, S) = 1 for all x. If we write the general solution to Eq. (|23p as 



y(x) = AC(x) + BS(x) 



(28) 
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then 

A = W(y,S), B = W(C,y) (29) 

This equation is quite useful for deriving relationships between the coefficients of the asymptotic expansions of the 
wavefunction in different regions of space as shown in sections ITT1 and HTT1 A more detailed discussion of the Wronskians 
is available in Powell and Crasemann's book on quantum mechanics^. 
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FIG. 1: Functions Ci(x) (solid line) and S^Oc) (dashed line) for e — 1 and the Gaussian barrier with vq = 2 (solid line) 




FIG. 3: Transmission probability for three Gaussian barriers 
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FIG. 4: Transmission probability for the Gaussian well as a function of the well depth for three values of the energy 




